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Abstract 

Single pion production off nucleons is studied in the framework of heavy baryon chiral 
perturbation theory to third order in the chiral expansion. Using total and some older 
differential cross section data to pin down the low-energy constants, most of the recent 
differential cross sections and angular correlation functions can be described as well as 
total cross sections at higher energies. We show that the contributions from the one loop 
graphs are essentially negligible and that the dominant terms at second and third order 
are related to pion-nucleon scattering graphs with one pion added. We also discuss the 
possibility of extracting the pion-pion S-wave scattering lengths from the unpolarized 
data. 
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1 Introduction 



1.1 General introductory remarks 

Single pion production off nucleons has been at the center of numerous experimental and the- 
oretical investigations for many years. One of the original motivations of these works was the 
observation that the elusive pion-pion threshold S-wave interaction could be deduced from 
the pion-pole graph contribution of the pion production reaction. A whole series of preci- 
sion experiments at PSI, Los Alamos, TRIUMF and CERN (and other laboratories) has been 
performed over the last decade and there is still on-going activity. On the theoretical side, 
chiral perturbation theory has emerged as a precision tool in low energy hadron physics. It not 
only allows to investigate low-energy pion-pion scattering to high accuracy, see e.g. refs.[|I], 
but also elastic pion-nucleon scattering and inelastic pion production can be studied in the 
threshold region, for a review see 0. Another important aspect of inelastic pion production is 
the excitation of resonances, some of which couple much stronger to the t^tcN final-state than 
to the pion-nucleon continuum. In what follows, we will be concerned with the low energy 
region which allows to address questions connected to the chiral structure of QCD. Since there 
already exist a few studies based upon various approximations, we first give an overview of the 
theoretical status and discuss in which respects these existing calculations can be improved. 



1.2 Short review of existing calculations and objectives 

Beringer considered the reaction t^N —>■ tttcN to lowest order in relativistic baryon chiral per- 



were 



turbation theory Low-energy theorems for the threshold amplitudes Vi and 
derived in These are free of unknown parameters and not sensitive to the vrvr-interaction 
beyond tree level. A direct comparison with the threshold data in the vr+p — > n'^TT'^n channel, 
which is only sensitive to Vi, leads to a very satisfactory description, whereas in case of the 
process TT~p — > vr^vr^n, which is only sensitive to V2, sizeable deviations were found from the 
total cross sections near threshold. These were originally attributed to the strong pionic final- 
state interactions in the J^^r = channel. However, this conjecture turned out to be incorrect 
when a complete higher order calculation of the threshold amplitudes T>i 2 was performed [Q. 
In that paper, the relation between the threshold amplitudes Vi and V2 and the tttt S-wave 
scattering lengths Cq and was investigated in the framework of heavy baryon chiral perturba- 
tion theory to second order in the pion mass (which is the only small parameter at threshold). 
The pion loop and pionic counterterm corrections only start contributing to the niiN threshold 
amplitudes at third order in the chiral expansion. One of these counterterms, proportional to 
the low-energy constant £3, eventually allows to measure the scalar quark condensate, i.e. the 
strength of the spontaneous chiral symmetry breaking in QCD. At that order, the largest con- 
tributions to Vi^2 stem from insertions of the dimension two chiral pion-nucleon Lagrangian, 
which is characterized by a few low-energy constants called q. In particular this is the case 
for the amplitude X>2- It appeared therefore natural to extend the same calculation above 
threshold and to compare to the large body of data for the various reaction channels. It was 
already shown by Beringer |^ that taking simply the relativistic Born terms does not suffice to 
describe the total cross section data for incoming pion energies up to 400 MeV in most chan- 
nels. In ref.0 the relativistic pion-nucleon Lagrangian including all dimension two operators 

^■^These are related to the more commonly used ^10 and A32 by A32 ~ \/TOI?i and Aio = — 2I?i — 3I?2, see 
also section 3.1. 
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was used, their coefficients being fixed from pion-nucleon scattering data and (sub)tliresliold 
parameters. Tliis parameter-free calculation lead to a satisfactory description of most existing 
total and differential cross section data as well as angular correlation functions. While being 
rather successful, that calculation could not give a definite answer to the questions concerning 
the importance of loop effects (and thus the sensitivity to the pion-pion interaction beyond 
tree level) and the convergence of the chiral expansion. Also, the amplitudes calculated from 
the tree graphs in ref.|^ are, of course, purely real. How severe this approximation is can 
only be judged after a one loop calculation above threshold has been performed. A ffist at- 
tempt to include third order contributions above threshold in the framework of heavy baryon 
chiral perturbation theory was described in ref. There, the pertinent one loop graphs and 
insertions from the dimension three Lagrangian in the formulation of ref.0 were calculated. 
However, only three of the five physical channels were considered and the finite pieces of the 
unknown dimension three low-energy constants were not determined but rather varied within 
given generous bounds. It was claimed that the contribution from the chiral dimension three 
amplitudes are large and play an essential role in the description of the data. Note that one 
of the low-energy constants (LECs) varied in ref. P] was shown to have a fixed coefficient in 
ref. |]I0[. Also, no explicit formulae for the various third order contributions were given, so that 
it is not possible to check the consistency of the amplitudes constructed there with the thresh- 
old amplitudes obtained in |Q. For these reasons, we do not consider the results obtained in 
ref. P] as conclusive. In addition, there exist some model calculations which are only partly 



constrained by chiral symmetry, see e.g. 12|. The one closest in spirit to a chiral expansion 



is the one of ref. [0, in which Beringer's Born terms where supplemented by explicit A and 
Roper (tree) contributions. Clearly, the inclusion of the resonances as done in that paper is 
not based on a consistent power counting scheme but rather it is argued that phenomenology 
demands the extension of the effective Lagrangian to include these higher mass states. 

Our objective is to perform a complete third order calculation using the minimal effective 
pion-nucleon Lagrangian. In addition, we will use all available information about the appearing 
LECs, in particular from the comprehensive study of elastic pion-nucleon scattering described 



in ref. [|T0|. Throughout, we will work in the so-called standard scenario of spontaneous chiral 
symmetry breaking, i.e. assuming a large scalar quark condensate. While the extension to the 
small condensate case is in principle straightforward, we do not consider it in what follows. In 
particular, our study allows to address the following questions: 

i) How sensitive are the data in the threshold region to the vrvr interaction beyond leading 
order? Previous investigations seem to indicate that a better method to access the vrvr 



interaction is the use of Chew-Low techniques (see e.g. [Q) rather than a direct use of 



the chirally expanded nnN threshold amplitudes (see e.g. refs.^, |T5|)- 

ii) How quickly does the chiral series converge? The previous analyses have not yet given a 
unique answer to this, although the study of ref. |p seems to indicate a slow convergence. 
On the other hand, the results obtained in the relativistic tree level calculation in ref. 
could be taken as an indication that loop effects are not very important. 

iii) At third order, new LECs appear. These can be determined from a ffi to cross section 
data. Are their values of natural size? If that is the case, it would also be interesting to 
develop an understanding of the numerical values as it was done for the dimension two 



LECs in ref.[|T6|. Furthermore, in ref.[§| resonance saturation was used to pin down the 
order three LEC contribution to the threshold amplitude D2 and this procedure can be 
made more precise when all LECs are determined from above threshold data. 
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iv) Can the existing data be described consistently in chiral perturbation theory? So far, 
mostly total cross section data and some angular distributions have been studied (with 
the exception of the detailed angular correlation functions in the ir'p Tr^vr+n channel). 
With the new TRIUMF data 



15 



on TT^p 



'n 



©HI, 



and the more recent ones from CHAOS 
this data base has considerably increased and allows for detailed tests 
on chiral pion-nucleon dynamics. In particular, we now have differential cross sections 
with respect to the invariant dipion mass squared M^^, the squared momentum transfer 
to the nucleon and the scattering angle between the two negative pions in the dipion 
restframe [T^. 



1.3 Organization of the paper 

This paper is organized as follows. In section 2 we briefly review the effective Lagrangian un- 
derlying the calculation. Section 3 contains some formal aspects including the definition of the 
T-matrix and of the pertinent observables. The chiral expansion of the invariant amplitudes is 
performed in section 4. Section 5 contains the results and discussions thereof. In the appendix, 
we give the explicit expressions for the amplitudes. 



2 Effective Lagrangian 

At low energies, the relevant degrees of freedom are hadrons, in particular the Goldstone bosons 
linked to the spontaneous chiral symmetry breaking (for a review, see e.g. We consider here 
the two flavor case and thus deal with the iso-triplet of pions, collected in the unitary matrix 
U{x) = M^(x). It is straightforward to build an effective Lagrangian to describe their interac- 
tions, called /^TTTT- This Lagrangian admits a dual expansion in small (external) momenta and 
quark (meson) masses as detailed below. Matter fields such as nucleons can also be included in 
the effective field theory based on the familiar notions of non-linearly realized chiral symmetry. 
The pertinent effective Lagrangian called Ct^n consists of terms with exactly one nucleon in the 
initial and the final state. The various terms contributing to a process under consideration are 
organized according to their chiral dimension, which counts the number of derivatives and/or 
meson mass insertions. Here, we work to third order in the corresponding small parameter q 
(which is a generic symbol for an external momentum or pion mass). Consequently, the effective 
Lagrangian consists of the following pieces: 

where the superscript [i) gives the chiral dimension. The relevant terms of the meson La- 
grangian read 

4^2 = ^{V^UV^U^ + x+) , x+ = u^XU^ + ux^u (2.2) 
4tl = |(V,f/V^f/t)2 + |(V,t/V.f/t)(V'^f/V^t/t) + ^(x+)^ 

+ ^{2(V^t/V^f/t)(x+) + 2(xWf/ + xf/W)-4(x^x)-(x+)'} + --- (2-3) 
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where the SU(2) matrix-valued field U{x) collects the iso-triplet pions, 



U{x) = 



\J F'^ — -K {xY + ir ■ TT (x) 



(2.4) 



with F the pion decay constant in the chiral limit and V^^ = (9^ + ... is the pion covariant 
derivative. Here, we only need the partial derivatives since we are not concerned with left- 
or right-handed external vector currents. This representation of the pion fields, the so-called 
cr-model gauge, is of particular convenience for calculations in the nucleon sector. The quantity 
X contains the light quark mass = = mp^ and external scalar and pseudoscalar sources. 



The ellipsis in eq. (|2.3| ) stands for other terms of order which do not contribute in our case. 
The fourth order terms in the meson Lagrangian serve to cancel some of the divergences of 
the loop diagrams and contain the mesonic low-energy constants ^1,2,3,4. The latter encode 
information about the chiral corrections to the vrvr scattering lengths. The finite pieces of 
the low-energy constants ii in eq.( |2.3| ) are renormalization scale dependent and are related to 
the ii of ref. [|18| via (since both of these sets are used in the literature, it is pertinent to give 
these relations) 

Ii = 16 tti TT^ (A) — 2 In — , ai = 6 , 0^2 = 3 , 0:3 = — 4 , 0^4 = 1 , (2.5) 

A 

and their actual values will be discussed later. Here, A is the scale of dimensional regularization. 
Let us now discuss in some more detail the terms appearing in the various parts of the pion- 
nucleon Lagrangian, i.e. in C)^^' . We make use of baryon chiral perturbation theory in the 
heavy mass formulation (jTU], ^ (HBCHPT). The nucleons are considered as extremely heavy. 



This allows to decompose the nucleon Dirac spinor into "large" (H) and "small" (h) components 

^(x) = e-'"^''"{i/(x) + h{x)} , (2.6) 

with Vn the nucleon four- velocity, f ^ = 1, and the velocity eigenfields are defined via i/ H = H 
and 1/ h = —h^^ Eliminating the "small" component field h (which generates 1/m corrections), 



the leading order chiral vrA^ Lagrangian reads 

C'^^l = H{ivD + gAS-u)H . (2.7) 

Here, = <9^ + L^ denotes the nucleon chiral covariant derivative, S*^ is a covariant gener- 
alization of the Pauli spin vector, the nucleon axial vector coupling constant (in the chiral 
limit)^ and = iu^'W ^JJu^ . To leading one-loop accuracy, i.e. order 0{q^), one has to 
consider tree graphs from 

£eff = + -Cfj^ + dfl; (2.8) 

(2) 

and one-loop graphs with dimension one insertions only. The pertinent terms of C)^{^ read 

+ (c2 - ^) (^ ■ + C3U-U+ (^C4 + ^) [S^, S'']u^u, + ...^H . (2.9) 



^^Throughout we will work in the isospin limit apart from kinematical corrections to be discussed later. 

#^Tlie role of is to single out a particular reference frame, here the tt°'N center-of-mass frame. 

*^We do not display explicitly the "°" commonly used to denote quantities in the chiral limit. The difference 
between the physical and the chiral limit values has to be kept in mind since it enters the renormalization 
discussed below. 
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12) 

All terms in are finite. The corresponding LECs have been determined in refs.[0, |21|, |10 



A physical interpretation of their numerical values can be found in ref. |T^, where it is shown 
that resonance saturation (mostly due to the A-resonance) can indeed explain the size of these 
LECs. The first divergences appear at 0{q^) in HBCHPT, the corresponding determinant has 
been worked out by Ecker p3|. We use here the minimal form of the dimension three irN 
Lagrangian as constructed in ref. ||T0[| 

23 _ 31 _ 



■-■kN 



i=l 



i=2A 



(2.10) 



The first set contains terms which stem from the 1/m expansion of the various dimension one 
and two operators and thus have fixed coefficients. Then there are 23 terms, from which 14 are 
divergent. For our calculation, we need the following dimension three operators 0\ 



(3) 



o. 



(3) 
ct 



1 + 8m C4 
32m2 



= idi{\)[u^,[vD,u^]]+i[d2{\) 

+ *(^3(A) + ^^ [vu,[vD,vu]]+i[d^{\)- 

+ d^{X) [x-, v-u] + dio(A) S-u {u-u) + dii{X) S'^u" {u^Uiy) 
4^A(l + 4mc4) + c/i^ 



9a \ 



disiX) + 



32m2 

2^^(l + 4mc4) 
IGm^ 



s-u{{vuy) 

S^v-u (u^v-u) 



+ c?i4(A)e'^^"V^^(b-/^,wjM,) + id,,+ 



9\ 



For these divergent operators, we introduce scale-independent renormalized LECs, 



di = di - 



L(A) + (47r)2lnA^ 



with 



L(A) 



ln(47r) + r'(l) + 1 



(2.11) 



(2.12) 



(2.13) 



167r2 U-4 

Note also that the operators Oi and O2 (0i4 and O15) appear as a sum (difference) for ttN — »• 
tcttN. Four combinations of counterterms {di + d2, d^, d^, ^14 — ^15) have already been fixed from 
pion-nucleon scattering data and dis from the so-called Goldberger-Treiman discrepancy. In 
addition, there are eight terms which are only necessary for the renormalization, and thus do 
not appear in matrix elements of physical processes. To be precise, these terms stem solely 
from the divergent part of the one-loop generating functional and have no (finite) counterparts 
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in the relativistic theory. They read 



0§l = d24{\)i{v-Df + d25i\)v-D S-uvD + d2ei\){i{u-u)v-D + h.c.) 
+ d27{X) {t {{v-uy) V-D + h.c.) + d28{X) {i vD + h.c.) 
+ MX) (S'^iv-D, u^]vD + h.c.) + MX) [e'^^'^^VaSpiu^, u,]vD + h.c.) . 

(2.14) 

All these terms are proportional to the nucleons equation of motion, ^ v ■ D H , and thus can 



be transformed away by appropriate field redefinitions. For more details, see ref. ||10 



3 Single pion production: Formal aspects 

In this section, we outline the basic technical framework for the reaction irN — > uttN above 
threshold treating the nucleons as very heavy fields. 

3.1 Invariant amplitudes 

We seek the T-matrix for the process 

7r'^(gi) + N{mv + pi) Ti\q2) + tt^I^s) + N{mv + pa) , (3.1) 

with N denoting a nucleon (neutron or proton) and vr*^ a pion of cartesian isospin a. Note that 
this process is characterized by five independent four-momenta (or five Mandelstam variables). 
It is most appropriate in terms of the chiral expansion to choose as variables the energies of the 
out-going pions and the three invariant momentum transfers (z = 1, 2, 3) defined below. Since 
we work in the framework of HBCHPT, the nucleon four-momenta are written in terms of a 
conserved four-velocity v and small residual momenta pi,2, where small means v ■pi^2 ^ ^ and 
m is the nucleon mass. The vrA^ —>■ txtxN transition matrix element can be expressed in terms 
of four invariant amplitudes, denoted A, B, C and D 0, which depend on the five momenta 
Pi,P2,qi,q2 and g3p3 as 



Tf,^ = {NAmv+p2)n\q2)n''iqs)\T\Ns{mv+p,)7T%q,)) 



M1U2XS 



S-qrA + S-q2B + S-q^C + le^^^^qW^q'^v" D 



abc 



Xs , (3.2) 



with Xs the conventional two-spinor and 



M^ = \P^. (^ = 1,2), (3.3) 
V 2m 

the standard spinor normalization factors which are mandated by the matching to the rela- 
tivistic theory if one works in the so-called Pauli interpretation of the heavy nucleon fields (a 



detailed discussion of this topic can be found in ref. ||2^). The invariant amplitudes are complex 



functions of the momentum variables since in some of the one-loop graphs considered here, the 



^®In what follows, we will suppress the dependence of the invariant amplitudes on the momenta. 
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intermediate pion-nucleon state can go on-shell. The pertinent tree graphs lead, of course, to 
real amplitudes. The isospin decomposition of the invariant amplitudes reads 



X 



abc 



X e {A,B,C,D}. 



(3.4) 



The following calculations can be simplified considerably if one notices that for most diagrams, 
the amplitudes have the following symmetries under exchange of particles (note that for these 
equations to hold, the momenta pi and qi have to be chosen in-coming and ^2, Pi and ps 
out-going) : 



^i(?i, 52,^3) 
^2(^1, ^2,^3) 



M{qi,q2,qz) = 
-Di(gi,g2,?3) = 



-52(-g2,-gi,g3) 
-C3(-g3,g2,-gi) 
^3(gi,g3,g2) 
-Si(-g2, -91,^3) 

-Cl(-?3,-9l,?2) 

C2(gi,g3,g2) 

-B4(-g2,-gi,?3) 

-C4(-?3,-gi,g2) 

-^2(-?2,-?l,?3) 

-Di{-qi,q2,-qi) , 



(3.5) 

(3.6) 
(3.7) 

(3.8) 

(3.9) 



and in all these expressions the fixed values of pi^2 are not made explicit. Thus for these 
diagrams it is enough to give the formulae of six amplitudes Ai, A2, A4, B3, Di and D4 in 
order to construct the full T-matrix element. However, in the tree graphs with a dimension 
two or three NNtth vertex, these symmetries are not present. In these cases, there appear 
expressions proportional to S ■ pi and S ■ p2- In the center-of-mass (cm) frame these equal 
—5" • gi and —5" • (52 + ^3) respectively, which obviously leads to an incquivalcnt role played by 
the pion momenta qi, q2 and ^3. For such diagrams the amplitudes Bi, B2, -84,-02 have to be 
given in addition to the six discussed before in order to fully determine the T-matrix element. 
The five experimentally accessible channels follow from the isospin amplitudes. 



up - 
7i~^p 
7i~p - 



71 p ^ 7r°7r°n 

7l~p TT^TT^p 



TT TT n 

> ti^tt'^P 



TT TT n 



X = V^(X2 + X3) , 

X = X3 + X4 , 
X = V2{Xi + X2) , 
X = v^Xi , 
X = X2 + X4 , 



(3.10) 
(3.11) 
(3.12) 
(3.13) 
(3.14) 



for X G {A, B, C, D}. Wc will calculate the amplitudes Xi in the isospin limit with the charged 
pion mass (M^ = 139.57 McV) and the proton mass (m = 938.27 Me V). Isospin breaking is 
accounted for in a minimal way by shifting the kinetic energy (in the laboratory system) of the 
incoming pion, called T^, from the isospin symmetric threshold 



'thr,iso 



mJi + 



3M^ 
'2m 



170.71 MeV 



(3.15) 



to the correct threshold of the corresponding reaction. For the five channels given above, the 
corresponding shift is ^T^r = +1.68, —5.95, +1.68, —10.21 and —5.95 MeV in the same order 
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as in eqs.( |3.10"| - |3.14|) . Finally, we need the expressions of the two threshold quantities T>i 2 in 
terms of the invariant amplitudes defined above. At threshold, ^2 = ^3 = and there are only 
two amplitudes. 



V2 



__^thr ^thr 



(3.16) 



with Af^^' = 1.006. The chiral expansion of Vi and has already been considered to third 
order, i.e. including all terms of order since the pion mass is the only dimensionful quantity 
at threshold and the lowest order terms are non- vanishing in the chiral limit. To end this 
paragraph, we note that from now on we use the energies of the outgoing pions and the three 
invariant momentum transfers squared as kinematical variables, because 002 = v ■ q2, uj^ = v ■ = 
0{q) and ^1,^2,^3 = C^(?^) have the expected chiral dimensions and consequently the residual 
nucleon energies v ■ pi are of second order. 



3.2 Observables 

To calculate the total (unpolarized) cross section, we need the invariant matrix-element squared 
multiplied by the appropriate weight functions. It reads 

s,s' 



- A* A {ujI - ql) + B*B {ul - qj) + C*C (cu| - ql) + {A* B + AB*) (cuiCUa - gi ■ gs) 
+ {A*C + AC*) (CU1CU3 - gi ■ ga) + {B*C + BC*) {uj2UJ; - q2 ■ q^) 

,2 „2\/„,2 „2\/„,2 „2\/i J2\f-\ I 1 {z — XiX2)'^ 



+AD*D {wf - qi){wi - qi){wi - qi){l - xi){l - xi) \l 



(l-x?)(l-xi)^ 



(3.17) 



where the angular variables Xi 2 and z are defined below. The formulae for total and the double 
and triple differential cross sections then become: 



0"tot(^7r) 

(Pa 
duj2dVL2 
d^a 
dcuodflodfl':) 



AmS 



(47r)5^|gi 
4m2|g2||g3|5 

(4^0V^¥S 



duo2duj^ / dxi 

2<1 J-l JO 



\M\' 



/ du^ / d 


i)\M\ 


Joj~ Jo 





\M\ 



(3.18) 

(3.19) 
(3.20) 



with s = (m + M,r)^ + 2mT.^ the total center-of-mass energy squared and Ui = 

is a Bose symmetry factor, 5 = 1/2 for identical pions and 5 = 1 otherwise. Here, (j) 
is the (auxiliary) angle between the planes spanned by g2 and gl as well as g2 and gjj. In 
accordance with the experimentalists convention, we have chosen the coordinate frame such 
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that the incoming pion momentum qi defines the z-direction, whereas q2 hes in the xz-plane. 
The polar angles 9i and 62 of the outgoing pions (with xi = cosOi) are in general non- vanishing 
and so is the azimuthal angle (p2 of n'^. By construction, the azimuthal angle of is zero. 
Furthermore, 



^3 = ^3 1 + 



^/s UJ2) + MI{UJ2 + U;3 - a/s) 



ul - Ml 



(3.21) 



The cosine of the angle between the two three-momenta of the outgoing pions, q2 and gs, 
respectively, is given by 

z = cos 9i cos 02 + sin 9i sin 62 cos (p2 , (3.22) 
and it can be used to express the energy as 



1 



Us 



2[{^ - UJ,y - z^\q2\ 



-|(v^-cj2)(s-2v^wi -m^ + 2M2) 



-z\q2\ 



2^/sUJl 



m 



2^2 



4M2(m2 + (1-^2)1^-21 



(3.23) 



In the formula for the total cross section, eq.( p.l8| ), the restriction 2;^ < 1 is equivalent to 
UJ2 G [IiJi], UJ3 G [u^,ujt], with 



It 



2{s - 2v/icu2 + M2) 



(Ti - M^f -m'^ + Ml 
2(v/i-M^) 

(v^ - 0J2){s - 24IUJ2 -m^ + 2Ml) 



±1^2 Iv (s - 2^/suj2 - m^Y - Am?Ml 



(3.24) 



(3.25) 



To compare our calculation to the most recent TRIUMF data, other combinations of differen- 
tial cross sections are needed. Consider first the differential cross section with respect to the 
invariant mass of the final dipion system. 



da 



mS 



dMl 



(47r)4^|gi \\q2 + | JT^{T^ + 2M^ 



dt / duj2 



(3.26) 



The expression for d'^a / {dMl^dt) can be easily deduced from eq. (|3.26| ). Here, the integration 
boundaries of t are given by 



t± = Ml + Ml^ - 2cuicu23 ± 2|gi| |g2 + g3| 
and Ci;23 is fully determined by M^^ via 

Ml, + s-m^ 



LV23 — UJ2 -\- u^ 



2^s 



(3.27) 



(3.28) 



Consequently, uj2 is restricted to the overlap of the intervals [Jj , 

r± "^23 



1, 2) with 



± 



(^i3-M2J(M2^-4M2) 
4M2 



(3.29) 
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Denoting by a the angle between gl and g2 + (h and by /5 the one between g2 and g2 + , the 
azimuthal angle 0' is defined via 



x\ = cos a cos /3 + y (1 — cos^ a)(l — cos^ /3) cos0' . (3.30) 
The differential cross section with respect to t = (gi — g2 — Qs)^ reads 
rfo- m>S r dMl 



The kinematically allowed region for M^^ is 



M^^ = A7<|t/2(s + m2-M2)+m2M2±|gi|^-st(4m2-t)[> , (3.32) 



and also M^^ e (-^i — m)^). Finally, the last measured quantity is the single differential 

cross section with respect to the scattering angle of the two outgoing pions, in the center-of-mass 
frame of the outgoing two-pion system: 

da 2mS 1 /-^^ . \M\^ 



dcose (47r) Vr^(r^ + 2M^) v^l^il J ""J J sin(^' 



=— -- / dM^ [ dcosa [ duj2 
AqA J J J 



Wl\W2\ 



X ^ (3.33) 

1^2! v (1 — cos^ a)(l — cos^ (3) 



with q'l = (gi)cM and ^2 = (^2)cm the pion momenta in the dipion cms, with magnitude 
= {Ml - 2Mlit + Ml) + {t- Mlf)/AMl^ and = (M^ - AMI)/ A respectively. 
The necessary kinematic constraints are given by the two cosine functions. This completes the 
necessary formalism for our calculation. 



4 Chiral expansion 

The effective Lagrangian described in secti on can now be used to work out the chiral expansion 
of the invariant amplitudes defined in eq. (|3.^ . Before discussing the explicit contributions, we 
first display the general structure of the chiral expansion. All the invariant amplitudes take the 
form 

^ _ ^tree ^ jS^loop ^ (4.1) 

where the tree terms are of dimension one, two and three whereas the loop graphs are of third 
order only. Corrections not calculated here start at order q^. The corresponding tree graphs 
are shown in figs. 1,2. Here, diagram 3p denotes a genuine third order counterterm specific to 
pion production and the graph 3s embodies the next-to-leading order pion-pion interaction. 
We note that many of the diagrams are simply elastic pion-nucleon scattering graphs with one 
additional pion attached in all possible topologies. The one-loop graphs are split into tadpoles 
(tl-tl4 in fig.3) and so-called self-energy diagrams (sl-s35 in figs. 4, 5). A typical one-loop vrvr 
interaction is e.g. tl3. We note again that many diagrams can be split into a vrA^ graph with 
an additional pion. This leads one to expect that a precise description of elastic pion-nucleon 
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scattering is an important ingredient to the calculation performed here. As an important check 
of our calculation we recover the explicit expressions for the two threshold amplitudes worked 
out to this order in ref. 0. The contribution from the tree graphs proportional to the dimension 
three LECs was estimated via resonance saturation in ref . . We are in the position to test this 
assumption (see below). Of course, many of the loop graphs shown in figs. 4, 5 simply amount 
to mass and coupling constant renormalization. This is by now a standard procedure and we 
refrain from discussing it here in detail (see the appendix). Also, the appearing divergences 
~ l/((i — 4), with d the number of space-time dimensions are of course cancelled by the 
appropriate infinite contributions from the dimension four mesonic LECs and the dimension 
three pion-nucleon LECs, first given in an overcomplete basis by Ecker p3], see section 2. This 



procedure closely parallels the calculations described in detail in ref.0. It serves of course as 
an important check of the calculation. 



5 Results and discussion 

We are now in the position to analyze the existing data within the framework laid out in the 
previous section. We first discuss the fits which determine the LECs and then show predictions 
and analyze in detail the chiral expansion, in terms of sensitivity to the tttt and vrA^ interactions, 
its convergence and related issues. 



5.1 The fitting procedure 



As stated in section 0, we have to deal in total with 19 independent combinations of LECs. 
take the values of the four mesonic LECs from ref. ||2^ , 



We 



^^1 GeV) = -5.95 , £^(1 GeV) = 4.35 , £^(1 GeV) = 1.64 , GeV) = 2.29 



(5.2) 



in units of 10~^. Although we did not perform a systematic calculation within the framework 
of generalized CHPT (i.e. assuming a small value of the quark condensate), we can get an idea 
about that approach by setting ~ —70 ||T8|, i.e. £3(1 GeV) = 117 ■ 10~^. Clearly, this is by 
no means a substitute for a complete calculation but should allow us to discuss some trends. 
A full analysis based on GCHPT is, however, not the topic of this paper. Turning now to the 
pion-nucleon sector, we note that the four dimension two LECs as well as five (combinations 
of) dimension three meson-baryon LECs can e.g. be taken from Fit 2 of ref.||10||. 



Cl 

di + J2 



-1.42 , C2 
3.31 , d-x = 



3.13 , C3 
-2.75 , J5 



-5.85 , C4 = 3.50 



-0.48 , dii-d 



'15 



-5.69 , di8 = -0.78 . (5.3) 



The values for the q and di are in GeV~^ and GeV~^, respectively. Therefore, we are left 
with six genuine LECs which have to be fitted. These are J4, dio, du, di2, di3 and die. ^oi 
the fitting procedure, we only use data with T^^ < 250 MeV. This rather large value is mainly 
motivated by the fact that in the very near threshold region, say in the first 30 MeV above 
the respective thresholds, there are essentially no data in the two channels 7r~p — > Ti^%~p 
and TT+p TT^Ti^p. The total cross section data are taken from refs.||25|-^^. Moreover, 



there exist double differential cross section data from Los Alamos [47, 48 1 and the Erlangen 
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group In ref.[^, a large amount of angular correlation functions are also given. 

All these data refer to the n^p vr+vr^n channel. In addition, there are the recent TRIUMF 
data listed in the introduction. We have performed a series of fits based on different input 
material. We refrain from describing all these in detail. What we will present here as results is 
based on fits to the total (including the recent TRIUMF data) as well as the double differential 
cross sections of refs. p7|, In that way, we can predict the angular correlation functions of 



ref. the new TRIUMF data |TJ] and the total cross sections for > 250 MeV. In these 
fits, we have found almost perfect anticorrelations between the values of the LECs diQ and di2 
as well as dg and dn. Therefore, we can only determine four independent LECs. In table |l], we 
give the results for various fits with all six LECs free and with one or two of them being fixed. 
We observe that due to the large anticorrelations, the fit with all six LECs left free leads to 



di 


None fixed 


dii = -5 


dii — diQ — — 5 


4 


0.26 ±2.29 


1.60 ±2.00 


0.97 ±2.00 


10 


-3.60 ±4.52 


-12.81 ±3.11 


fixed 


11 


-20.58 ±5.90 


fixed 


fixed 


12 


4.11 ±4.32 


12.25 ±3.27 


4.28 ±0.45 


13 


22.00 ±5.78 


6.60 ±0.36 


6.08 ±0.32 


16 


-4.70 ± 0.84 


-3.67 ±0.74 


-3.84 ±0.70 


xVdof 


2.20 


2.25 


2.29 



Table 1: The dimension three LECs di in GeV ^ from the various fits as described in the text. 



rather large values for some of the LECs. If one fixes two to any value of natural size, which 
in our normalization is of order one, the remaining LECs come out to be of natural size too, 
compare the last column in table ^ The uncertainties given in that table are the parabolic 
errors of the MINUIT package and have to be taken with some caution. In what follows, we 
will always use the values for the di as given in the first column of table |l]. In fig.p we show the 
fit to the total cross sections in all five physical channels and in fig.^ the corresponding double 
differential cross sections. We note that there are some inconsistencies between the older and 
more recent data, most pronounced in the two channels Ti^p Tr+vr+n and n^p vr+Tr^n. 
Note also that the data of the OMICRON collaboration ||33|, |2^ have been criticized concerning 
the normalization. If we perform fits without these data, the x^/dof only improves by one 
permille. The LECs di change mildly, the only appreciable difference is found for d^. This LEC 
only appears in the amplitude D^, which contributes to the channels with one neutral pion in 
the final state. The data from the OMICRON collaboration feature prominently in the process 
7r~p — > Ti^Ti~p which explains the sensitivity of the LEC d/^. In what follows, we keep the data 
of ref. 26 1 in our data base. We have also performed fits using the LECs q and di from fits 1 
and 3 of ref.[jlO|. The resulting dimension three LECs only change moderately and the x^/dof 



is 2.21, 2.20 and 2.20 for fits 1,2 and 3, in order. None of the conclusions drawn in the following 
paragraphs depends on this choice, we always use the values from fit 2 given in eq.( ^.3|) . 



"^^As described in |^], neither chiral perturbation theory nor any of the resonance models can describe the 
double differential cross section data at = 1.301 GeV from ref.[^. We therefore decided not to use these in 
the fit. Also, we do not use the unpublished data of refs.pft pl[. 
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5.2 Predictions, further results and discussion 

We now turn to the predictions. In fig.§ we show the total cross sections for incoming pion 
momenta up to 400 MeV. Despite the large momentum transfers involved, the chiral description 
works fairly well, in particular in the 7r~p ir^n'n channel. On the other hand, for the two 
reactions n^p —>■ TT^n^n and TT~p —>■ n^Ti^p the chiral prediction overshoots the data for 



energies larger than 300 MeV. The angular correlation functions as defined and given in ref . [HQ 



are mostly well described, see figs.^,|TU|, with the exception of the smaller 6*2 values, a trend 
already observed in ref. 0. If one insists that these data are also fitted, the x^/dof worsens 
considerably. The comparison to the recent TRIUMF data is shown in fig.0 for da/dM^^ 
and da/dt. Whereas the data in the 7r~p channel are mostly well reproduced, there are some 
deviations for the vr+p process, most notably in da/dM^^. The description of d^a / dtdM^^ and 
da /d cos 6 is less good, we refrain from showing these. Again, if one includes these data in the 
fit, the x^/dof becomes intolerably large. This means that the presently existing data base 
shows some inconsistencies, a fact which also limits the precision of our description based on 
the chiral symmetry of QCD. 

We now consider the chiral expansion in more detail. We remark that in three channels the 
series seems to converge well {ir^p 7[^n~^n, Tr^p — >■ n^ir^p and 7r~p — > 7r°7r~p), in the two 
others (vr~p —>■ -k^-k^ti and Ti^p n^TT^n) the third order corrections become large even close 



to threshold. As two representatives of these classes, we show in the upper part of fig.|T2 
the close to threshold region (T^r < 210 MeV) for 7r~^p ix^-K^n (good convergence) and 
TT~p —>■ n^ir'n (poor convergence). However, it is of importance to further analyze these third 
order contributions. As shown in the lower part of fig.^ in all cases the contributions from the 
loops and (in most cases) from the terms proportional to the dimension three LECs are small 
(this holds also for the channels not shown in the figure). If the third order correction is large, 
it comes entirely from the 1/m corrections to the tree graphs proportional to the dimension two 
LECs Ci and the ones with fixed coefficients. The smallness of the unitarity corrections is a very 
important result. At first, it appears to be surprising since in elastic pion-nucleon scattering, 
the loop contributions are sizeable (in some of the invariant amplitudes) already in the threshold 
region. We note, however, that such terms appear with an additional pion line attached before 
or after the irN ttN subprocess, which leads to sizeable cancelations in the imaginary parts 
for the reactions considered here. This finding explains why the relativistic calculation based 
on dimension one and two tree graphs of ref.0 worked so well, since in that calculation all 
imaginary parts were neglected. Similarly, the commonly used resonance models like e.g. the 



one of the Valencia group |]TT| only acquire small imaginary parts from the finite width of 
the resonances, which is built in by modifying the corresponding propagators. Our analysis 
for the first time gives a reason why such a seemingly ad hoc procedure can work. Another 
consequence concerns polarization observables, which usually stem from the interference of real 
and imaginary parts of certain invariant amplitudes. In order to differentiate between different 
approaches to a given reaction, one often has to analyze polarization data, one prominent 
example being the multipole separation in pion photoproduction off nucleons. Due to the 
intrinsic smallness of the imaginary parts for inelastic pion production, one should find out in 
which polarization observables the small imaginary part is most efficiently amplified by the real 
part of a large amplitude. In principle, the framework outlined here can be used to do that, 
but due to the absence of polarization data, we refrain from discussing such an analysis here. 

We now consider briefiy the threshold amplitudes I^i,2, as defined in eq.( |3.16|) . In ref. |^ it was 



found that no tree graphs with dimension three LECs generated from resonance excitation (via 
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the dominant Roper resonance) contribute to "Di, whereas two distinct terms were considered for 
1^2 • The first is hnked to the deviation from the Goldberger-Treiman relation and corresponds 
to our term ~ rfig, the second one was the contribution due to the Roper decay A^*(1440) — *• 
N{7i7i)s, where the two pions are in a relative S-wave. The uncertainty related to this decay 
in fact was the largest source of theoretical uncertainty for the chiral description of V2. In our 
more general approach, we also have a contribution from the dimension three LECs to T>i. If 
we convert our central values and errors of the corresponding di {i 7^ 18) into the contribution 
to "^1,2, we find 

= (0.27 ± 0.07) fm^ , V^' = (-1.60 ± 0.15) fm^ , (5.4) 

which is within 1.5 a of the result of ref. 0] based on Roper excitation, X>^* = (— 0.40±0.90) fm^. 
In addition, we have a smaller contribution to Vi, which is however of the same size as the 
other contributions form the loops, q and the ii terms, cf. table 2 of ref. The theoretical 
uncertainties in eg. (|5.4|) have been obtained from the MINUIT uncertainties in the di taking 



into account the correlation matrix. 



Finally, we consider the sensitivity to the LEG £3. In fig.|T3| we show total cross sections for 
the two channels ir'^p TT'^ir'^n and ir^p TT^ir^n at T^^ < 220 MeV (which have already been 
investigated in the study of the low-energy theorems 0). In both cases the error bars of the 
recent data are smaller than the difference in the curves with £3 = 2.9 (standard case) and 
—70 (generalized case), in order. A similar trend is observed for the differential cross sections 
da/dt (at T^r = 220, 240 MeV) , which have been used in ref.|]n[ to extract by Chew-Low 
techniques. Our analysis tends to support the finding of ref.]^, namely that one can extract Oq 
from the threshold data, but with an uncertainty which includes both the standard as well as 
the generalized scenario. Presumably, a study of polarization observables would give a better 
handle on this question. In ref. 0, a smaller theoretical uncertainty on Oq was given, based on 
the observation that there is no A^*(1440) contribution to the threshold amplitude Vi. Taking 
into account the additional contribution given in eq. ( |5.4|) , we conclude that the theoretical 
uncertainty for Oq was underestimated in ref.p. 



6 Summary 

We have performed a complete third order calculation of the reaction vrA^ ttttN in heavy 



baryon chiral perturbation theory based on the minimal Lagrangian developed in ref. |]I0|- We 
had to consider 26 different tree graphs and 49 one loop topologies as displayed in figs. 1-5. 
Note that the tree contributions with fixed coefficients are obtained from the 1/m expansion 
of the relativistic amplitudes. The pertinent results of this investigation can be summarized as 
follows: 

i) We have used the total cross sections in all five physical channels (for T^- < 250 MeV, 
with the kinetic energy of the incoming pion in the laboratory frame) and the older 
double differential cross sections d'^a/dQdT for the process 7r~p — > 7j'~^7r~n to fit the six 
new dimension three LECs. The other four dimension two and five (combinations of) 
dimension three LECs were taken from the study of elastic vrA^ scattering in ref. |T^. We 
observe that the values of two pairs of these LECs are almost perfectly anticorrelated, so 
that only four LECs can be determined. They come out of natural size. 
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ii) Using this input, we predict the total cross sections for energies up to = 400 MeV. We 
find an excellent description of the ir^p tt^tt'ti channel whereas the largest deviations 
are seen for the process n^p vr+Tr+n, see fig.§. The angular correlation functions for 
TT~p — > 7r+'7r~n can be satisfactorily reproduced, with the exception of the small 62 angles, 
cf. figs.p,p!0|. In addition, most of the recent TRIUMF data on da/dM^^, da/dt, da /d cos 6 
and d^a/dtdM^^ can also be reproduced. 

iii) At third order, the contribution from the loop graphs is essentially negligible. Unitarity 
corrections therefore play no role. This allows one to understand why resonance models 
like in refs.[^, |T2[ work fairly well even in the threshold region (although these are not 
as precise as the calculation presented here). The effect of the terms proportional to the 
dimension three LECs is somewhat more pronounced. By far the largest contribution 
at this order comes from the 1/m corrections to the dimension two LECs q and the 

corrections with fixed coefficients. This together with the smallness of the loop 
contributions explains why the tree calculation in ref.0 works so well. It also means that 
by far the most important terms are the pion-nucleon subgraphs with an additional pion 
added in all possible topologies. 

iv) Our study is based on the standard scenario of chiral symmetry breaking where a quark 
mass insertion is counted as second order in the chiral expansion. This refiects itself in 
the natural size of the mesonic LECs ^1.2,3,4 used here. The generalized scenario with 
rUq ~ 0{q) can be modeled by setting £3 ~ —70 HTSj . p^^ Keeping all other mesonic and 



baryonic LECs fixed, we have studied the sensitivity of the total and differential cross 
sections as well as the angular correlation functions to the value of £3. We conclude that 
one is not able to pin down the LEC £3 with sufficient accuracy to discriminate between 
the two scenarios of chiral symmetry breaking by just comparing the observables directly 
with chiral analysis (which is of course different from applying Chew-Low techniques). 
We have also evaluated the dimension three counterterm contributions to the threshold 
amplitudes and discussed the resonance saturation approximation used in ref.p. 



There are two directions in which the study presented here should be extended. First, it 
would be interesting to study the sensitivity to polarization. In particular, it is conceivable 
that some polarization observables will be more sensitive to the pion-pion interaction than 
the unpolarized observables studied here. Second, a complete one-loop calculation requires to 
include the fourth order terms. For doing that, one first has to work out elastic pion-nucleon 
scattering to that order. Such investigations are underway and we hope to report on their 
results soon. 



*^°Clearly, this can not substitute for a detailed study within generahzed CHPT. Such a procedure should, 
however, be sufficient for simply getting an idea about the sensitivity to the tttt interaction. 
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A Contributions to the invariant amplitudes 



In this appendix, we explicitly give the contributions from the tree and loop graphs shown in 
figs 1-5 to the invariant amplitudes A, B, C and D. Only the minimal number of non-vanishing 
amplitudes as explained in section 3 is given. 

The calculation of the Born amplitude with fixed coefficients has been done by 1/m expansion 
of the relativistic amplitude. In this expansion, uj^ and ujz have been treated as quantities of 
order g, ti = (^i - g2)^^2 = (gi - gs)^ and = (g2 + gs)^ = count as 0{q^). As a 



consequence. 



V-P2 



{UJ2 + cosf - Ml 
2m 

2m 

(cU2 + 003) + 0{q^ 



+ ... 



+ 



(A.l) 

(A.2) 
(A.3) 



where the ellipsis denotes higher order terms. In diagrams 3a, 3b, 3n and 3o, there appear 
expressions va. S ■ pi = —S ■ qi and S ■ p2 = —S ■ (52 + gs)- As stated before, for these diagrams 
the amplitudes Bi,B2,B4,D2 have to be given in addition to the six usual ones, in order to fully 
determine the T-matrix element. 

1. Counterterm amplitudes 
Diagrams 2a+2b: 



Ai 



A2 = 
A, = 
Bs = 

Diagrams 3a+3b: 



■9a , . X 
-z-^v • (Pi +P2) 



~{-4ciM'^ - 2c2ijJ2^3 - 2c3g2 • ga) 



^c^{u2^3 - g2 ■ gs) ( -9 + -9 

■9A , .V-{pi+p2) 
1-=;:tC4{uJ2UJ3 - q2 ■ gs) o 



.9a , . 

%j;^Ca[uJ2U}3 - q2 ■ gs) 



2 2 , X / 1 1 

■ [pi-p2)\ — 2 

UJ2 LO3 \UJ^ LUi^ 



■ 9A f .V-{pi+P2) 

i-^Ci[uJiUJ3 - qi ■ qs) ^ 



. 9 A 



2F3 

''2F^' 



2 2 . J I I 
h v • (pi - P2) — 2 

(jj2 L03 W2 '^i, 



■9a , , / 1 1 1 



(A.4) 



^1 



A,. = 



—I 



■ 9 a 



mF^ 
: 9a 

mF^ 



C4[llIi{llI2 + 1^3) - coi - col - qi ■ (g2 + ga) + 2M^] 

[-4ciM^ + 2c2u;iu;3 + 2c3qi ■ qs - C4{u;iu;2 - 0^2 - qi ■ q2 + M'^)] 
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^4 = i^|^C4[a;i(a;3 - 0J2) + col - UJ^ - Qi ■ {qs - ^2)] 

Bi = i-^^[-4ciM^ - 2C2LU2CO3 - 2csq2 ■ qs - Ci{uJiUJ2 - - gi • ^2 + M'^)\ 

B2 = -^:J^C4[u;2(t^3 - t^i) + + - 52 ■ (gs - qi) - 2M^] 

-B3 = i^^^[-4ciM + 2C2OJ1U2 + 2c3gi • q2 + C4(a;2a;3 + - ^2 • ga - )J 

Bi i^|^C4[(a;i + a;3)(a;i + a;2 + a;3) - (gi + ga) • (gi + 52 + ga)] 

£'2 = i-^Ci 

(A.5) 

Diagrams 3c+3d: 



Ai — —I 



{ui +UJ2 + uj-iY - (gi + g2 + ga) 



2 



J{uJi+ uj2f + ujI - (gi + ga)^ - 
-Ci[uj20J2, - q2 ■ ga) 



-(-4ciM^ - 2c2a;2a;3 - 2c3g2 • g3) 



/ ,2 
LO2 



^ {ui + uj:if + ul - (gi + gg)^ - M^ ' 



a;. 



. . . (a;i+u;3)'+^i-(9i + ?a)'-M2 
^2 = -«7r^^4(w2Wa - 52 • ga) 2 

. ._9A , ^2)' - col - (gi + g2)' + 
^ -^^^4-^2C.a - ^2 • ga) -2 

{ui + - ul - (gi + g3)2 + M2\ 



D . 9a , + a;2 + u-i)^ - (gi + g2 + ga)^ 



n - 

-Di = -I- — PPTC4 



4mF3 

+ a;2)^ - - (gi + g2)^ + _ + 6^3)^ - o^l - (gi + ga)^ + 



n ■ 9a / + a;2 + a;3)^ - (gi + g2 + g3)^ 

— ^"i ?^'^4 



{ui + a;2)^ + - (gi + g2)^ - + u^f + ujI- (gi + g3)^ - 

"T r, + 



(A.6) 

Diagrams 3e+3f+3g: 

A = -^^(t^2^^a - g2 • ga) 
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1 



1 



{UM^dieiX) - GM^dis - J25(A)a;2C^3 + d29{X){ujl + uj^ + col)) 



+ 



UJ1UJ2 



(12M^di6(A) - 6l\rdi8 + d25{X)uJiUJ2 + o?29(A)(a;{ + 0^2+ uj^)) 



-i^A{uj2UJ3 - q2 ■ qs) 



2F3' 
1 



+ 



UJ1UJ3 
1 



(12M^die{X) - 6M^di8 + d25{X)u;iu;s + d29{X){u;l + ui + uj)) 



1 



+ 



(12M^c?i6(A) - 6M'di8 - d25{X)uJ2UJ3 + d29{X){u;i + + u^)) 
{12M^diQ{X) - 6M^di8 + d25{X)LJii^2 + d29{X){Ljl + + ujD) 



UJIUJ2 

iS^{uJiuJ-i - qi ■ qs) 



+ 



+ 



2F3 
1 

UJ1UJ3 
1 

UJ2OJ3 
1 



Da = 



0UiU}2 

.91 



{12M^diG{X) - 6M^di8 + d25{X)uJiUJ3 + d29{X){ujl + ujI + cuD) 



{12M'^die{X) - 6M^di8 - d2^{X)u2U3 + d29{X){ul + o;^ + uD) 
{12M^diQ{X) - 6M^di8 + d25{X)LJiUJ2 + d29W{i^f + + i^i)) 



4^3 

1 



1 



UJ2UJ3 

1 



(UM^dieiX) - 6M^di8 + d25(X)u;iu;3 + d29(X)(u;l + 00^ + uD) 
{12M'^die{X) - GM'^dis - d25{X)uj2UJ3 + d29{X){ujl + 00^ + cuD) 



+ 

UJl(jJ2 

A4 = Di = 
Diagrams 3h+3i: 



(12M2di6(A) - QM^di8 + d2^{X)uJiUJ2 + d29{>^){^l + + ^3)) 



(A.7) 



^3 

-^^(^2C^3 - ?2 • 53) 
1 



a;ia;3 



■((i24(A)(u;^ + + im^d28{X)) 



1 



■(ci24(A)K' + u;2') + 16M2d28(A)) ((i24(A)(u;2' + ^3) + l^M''d28{X)) 

UJ1UJ2 CJ2W3 



1 



a;ia;3 



-(d24(A)K' + + 16M2d28(A)) 



^3 

-«^(^2t^3 - ?2 • ?3) 

+ ^(ci24(A)K2 ^ ^2^) ^ 16M2J28(A)) + ^(J24(A)(u;2' + ^3) + 16M2J28(A)) 
a;ia;2 UJ2UJ3 
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A, 







-B3 = ^2^^^^^^ - ?i ■ ^3) 
1 



-{d2A{X){i^l+Ujl) + 16M'd28(A)) 



a;ia;3 

,2\ I ic7i/r2. 



+ (c?24(A)K + ^2) + 16M'ci28(A)) + {d2i{X){oji + ^3) + 16M"o?28(A)) 



Da 



—I 



:9\ 



1 



UJ1UJ2, 



-{d24{X){ool + cul) + l6M^d28{X)) 



+ {d24{X){0jt + OJ2) + 16M'rf28(A)) {d24{X){uji + iui) + 16M'(i28(A)) 

UJ1UJ2 UJ2i^3 



Di = 
Diagrams 3j+3k: 



(A.8) 



_^^^^^?_^(4M2(ii6(A) - 2M'dis] 



—I 



2F3 uji 
^ '^'^'^'i4M'd,,iX)-2M'd,s) 



2F3 UJ2 

A4 = Di = D4 = Q 



(A.9) 



Diagrams 31+3m: 



Diagrams 3n+3o: 



Ao = 



^3 

Ai 



_,lA^^l^^d2,{XH - 8M%siX)) 

^,lA.^^l±^^d2,{X)ul - SM%s{X)) 
Zr ^ UI2 



(A.IO) 



^1 = 'T^ 



gC2 uj2qi ■ 93 + ^3^1 • g2 _ ^ / cj3(cj| - M^) _^ a;2(a;| - M^) 



a;2 



<^3 



+8J26(A)g2 • qs + 8(i27(A)a;2<^3 + 8M^J28(A) - SdsoiX) (002003 - ^2 • ^3) 



^2 = 



C4 / {1^2 - (^3) (?2 • ^3 + M"^) + ^3^1 • ^2 - '^2^1 • ^3 , '^2{'^'i " 



m 



+ 



0J3 



+^di{X)+d2) 



{u;2 - LV3)q2 ■ q3 



+ MsiX) 



{U2 - LJ3)U2U}3 



8M^d5{X) 



OOi OOi 

+2{dii[X) - di5) h -d24{X)uJi{uj2 - UJ3) 

LO2 z 



-4M2J28(A) 



<JJ2 — 00-i 



+ Ad3o{X){LJ2l^3 - q2 ■ ^3) 



20 



D4 



- 9A 



i^zqi ■ q2 - i^2qi ■ qs , {^1 + cu3)(a;2t^3 - 92 ■ qs) + i^si^i - M^) 



+ 



0J2 



{Ui + lxl2){uJ2UJ2, - q2 ■ qs) + 1^2 (l^i " M"^) 



Bi — I— 



003 

C4 / -{uji + a;3) (gi • q^ - M^) - uj^qi ■ q2 + ujiq2 ■ qz ^ - M^) ' 

ml ijJ2 0^3 



F3 

+4(di(A) + d2) 



+a;3)gi • gs (cui + a;3)a;ia;3 o s^2 . /xn ('^i +^3) 



+2(di4(A) - dis) 



+ 44(A) - 
cc;2 a;2 

{u2 - UJ^) {uJiUz -qi-qz) 1 j 



+ 8M'4(A)- 



0J2 



- 2^24(A)a;2(a;i +a;3) 



+4M'J28(A) 



oJi + a;3 

UJ2 



+ 4J3o(A)(a;ia;3 - gi • q^) 



Bo = 



^C2 ujiq2 ■ q^ + i^zqi ■ g2 _<^( - M"^) ^ - M'^Y 

m i02 m\ijj\ Lij'i 



+8d26(A)gi • qz + 8d27(A)a;ia;3 - 8M^d28(A) - 8ci3o(A)(a;ia;3 - qx ■ qs) 



B. = 



.gA 



C4 / uJsiujf - M^) {uji + uj3){qi ■ qs - M^) + ujaqi ■ q2 - Uiq2 ■ qs' 



m 



+ 



OJ2 



-4(di(A) + d2) 



UJ2 



-4rf3(A)- 



<^2 



-8M'rf5(A)- 



c<;2 



-2(di4(A) - dis) h -a24(A)a;2(t^i +^^3) 



-4M'd28(A) 



+ Go's 

a;2 



+ 4j3o(A)(a;ia;3 - q^ • qs) 



{ix}2 + Uz){uJiUJz - qi ■ qs) + ^si^f - M"^) a;3gi • q2 - uJiq2 ■ qs 



+ 



UJ2 



{lui + LLi2){i^iUJ3 - qi • qs) + <^i('^i - M'^) 



■ 9A 
'2mF3'^^ 

. 9a 
'2mF^''^ 



UJ3 

LUi +UJ3 tUi+ UJ2 



<^2 



<^3 



'^2 + ^^^3 _^ ^1 + ^2 



<^3 



-6i^J3o(A) 



(A.ll) 



Diagram 3p: 



Ai 
A2 



i 

F3 



2—gAU2U3 - Adio{X)q2 ■ qs - 4(ii2(A)u;2C^3 - AM^dieiX) + 2M^di8 
m 

5'A<^2<^3 - 2dii(A)g2 • qs - 2di3{X)u;2UJ3 + 2M^cii6(A) - M'^dis 
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5, 



F3 



C4 



gAUJiUJ3 - 2du{X)qi ■ qz - 2cii3(A)a;ia;3 - 2M die (A) + M dis 



m 



L»4 = 

= 

Diagram 3q: 



+ ^f^29(A)(Wl - 

^12c^4(A) 
Di = 



(A.12) 



A2 
A, 



{qi - ?2 - - 

~F^(gi-g2-g3)'-M2 

i 1 

^(gi- 92 -93)^-^2 
= = 



[4M2rfi6(A)-2M2di8](t3-M2) 
[4M2di6(A) - 2M2di8](t2 - M^) 



[4M2(ii6(A) - 2M^dis]{ti - M^) 



Diagram 3r: Combine with diagrams Rle and s35. 
Diagram 3s: 



(A.13) 



A, = 



.9 a 



A, 



+4^2[g3 ■ (92 + 93 - 91)91 -q^ + qi- (92 + 93 - 91)91 • 93] 

+2M^£4[qi • (92 + 93 - 91) + 92 • 93 + (92 + 93)^ 
+4^2(93 • (92 + 93 - 91)91 • 92 + 9i • (92 + 93 - 91)92 • 93] 



-2M^£4[q2 ■ (92 + 93 - 9i) + 9i ■ 93 - (91 - 93)' 

Bs = _ _ [8^193 • (92 + 93 - 9i)9i • 92 

+44I92 • (92 + 93 - 91)91 • 93 + 9i • (92 + 93 - 91)92 • 93] 
-2M^e4[q3 ■ (92 + 93 - 91) + 91 • 92 - (91 - 92)^] 

A4 = Di^D^^O 



(A.14) 



2. Tadpoles 
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We use the loop functions as defined in ref.pT 
Diagrams tl+t2+t3: 



Ai = - — {uJ2UJ3-q2-q3)^r, 



Diagrams t4+t5: 



A2 = — - — {uJ2UJ-i-q2-q3)A„' 



4 \UJ1UJ2 UJ1UJ3 002^3 

32 , X A / 1 1 1 

B3 = — — (u;iu;3 - qi ■ gsjA^ 1 h 



4 \ci;iCi;2 ^1^3 ^2^3 



8 F^ ^ \uJ1UJ2 UJ1UJ3 
= Di = 



Diagrams t6+t7: 



4 i' ^ 

Ai = A4 = Di = D4 = i 



8 F5 " cui 



Si u;i + cug 

8 r° UJ2 
Ai = Ai = Di = D^ = 

Diagrams t8+t9: 



Diagram tlO: 



.gAUJ2 - UJ3 

F^ oji 
.gA UJi+uj3 

^ F5 002 

Ai = A4 = Di = D4 = 



A, 




2(t3 - M')Io{h) 

3gi-g3-| + 2M^)/o(t2)-4A^ 



167r2 



2M^ 
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5, 



1 

+3 

i_9A_ 

1 

~3 
4F5\3 



(3q, • ?2 - I + loih) - 4A. + (sM^ - |)] | 

{2M^ - h)Io{t3) - 2K + (3^2 - 2M^)Io{t2) - 2A^ 



3^2 • ^3 - I + 2M^^ loits) + 2A, + ^ 



2M' 



^3 



_(2m2 - ^)/o(t2) + 2A, - (3ti - 2M^)Io{h) + 2A, 



Diagram til: 



(A.19) 



Diagram tl2: 



Diagram tl3: 



A4 



A2 
A4 



9A 



4 F^ ^ 
8 " 

2li gA^ 

8 



2F5 (gi - 52 - ga)' - M2 



—I 



1 



2F5 (gi - g2 - gs)' - M2 



A,(i3 - M^) 
A,(t2 - M') 



2F5 (gi - 52 - ga)' - 
= D4 = 



(A.20) 



(A.21) 



Ai 



'2F5 (5, - 52 - 53)' - 
|(t3 - M^)[Io{t3){2ts + M^- 2qi ■ (52 + ^3)) - 4A,] 
+7o(i2)[M2(2M2 - ^2) + ^^(4^2 - lOM^) - ^^(2^2 - 2M^)] 

+A,(-4M2 + ^q, . 52 - ^q2 ■ qs) - 7^(6^' " ^2)^2 • (^i + ^3) 

+Io{ti)[M\2M^ - + ^^{iti - lOM^) - ^^(2^1 - 2M^)] 

3 3 
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A, 



1 



2F5 (gi - g2 - gs)^ - M2 



/o(t3)[M2(2M2-t3) + 



^(4t3 - lOM^) + ^(2t3 - 2M2)] 



+A,(-4M2 + -gi • g2 + -gi ■ ga) - ^ 
+ {t2 - M^)[h{t2){2t2 + - 2g2 ■ (gi - gg)) - 

'?2 ■ gs 



(6M2 - t3)gi ■ (g2 - ga) 



. ^ H2- h;^^^^ _ ^Q^2) ^ ^i:^(2ti - 2M2)] 



+Jo(ti)[M'(2M'-ti) 



+A^(-4M2 - -g2 ■ ga + -gi ■ ga) + J^""^' 



(6M2-ti)ga-(gi + g2) 



9a 



1 



2F5 (gi - g2 - g3)2 - 



Ioih)[M\2M^~t,) + 



+A^(-4M2 + ^gi . ga + ^gi ■ g2) - ^(SM^ - ta)gi ■ (g3 - g2) 



r2,o..2 .^ ?2-_g3^4^^_^Q^2)^gi-g2,o. 0.^2, 



+/o(t2)[M^(2M'-t2) - 
+A,(-4M2 - ^g2 ■ g3 + ^gi ■ gs) + 



^(2t2-2M2)] 



(6M2 - t2)g2 ■ (gi + gs) 



A, = 
Diagram tl4: 



+ {ti - AP)[Io{h){2t, + - 2g3 ■ (gi - gs)) - 4A^] 



(A.22) 



Ai 

A2 

5a 
A, 



■ 9 a 



1 



>Mgi - g2 - ga)' - " 

1 

(gi - g2 - ga)' - 



A, 
A, 



"F^ (gi-g2-g3)'-M2^'^ 

= L)^ = 



-5ta + ^M' - (g2 + ga - gi)' 
-5t2 + ^M' - (g2 + ga - qif 
-5ti + ^M' - (g2 + ga - Qif 



(A.23) 



3. Self-energies 



We use the loop functions as defined in ref.pl 
Diagrams sl+s3: 



Ai = (^y-^ {UJ2UJ3 - q2 



- g2 ■ gs) (^5' + 



J2M - J2(-CUl) (]_^]_ 
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J2(a;2) - J2{-i^2) / 1 1 



, ,2 
<^2 



) + 



to- 









0J2) _ 



J2(a;i) - J2(-a;i) _ 

.U2 U2,, 



J2M - M-^^) fj_ J_\ J2(C^3) - J2(-^3) 



<^2 



a;3 



1 1 

h — 

OUi UJ2 



A, 



J2M + J2{-0J2) - 2J2(0) / 1 1 



J2M + J2(-C^l) - 2J2(0) 



/ ,2 
^1 



1 

0^2 



J2(C^3) + J2(-a;3) - 2J2(0) 



1 ^ 1 



-^2 ((^2) - J2{-^2) /II 



^1 



^ i2 
<^2 



+ 



^ ,2 



/ ,2 
^1 





1 


VCJ2 


h — 


1 

h 


a;2/ 


1 


) 


<^3 





J2(a;3) + J2{-0Jz) - 2J2(0) 



^ J2(C^2) + J2(-C^2) - 2J2(0) /J_ _ l_\ 



, ,2 
^3 



8 vf; 



^2 ^3 / 



^ J2('^2) - ^2(-C^2) /J \_\ ^ J2{i03) - ■/2(-^3) 



1 1 
a;2 



Diagram s2: 

A2 
A, 



(y) (a;2a;3 - q2 ■ Qs) (s^ + ^) 



2a;ia;2a;3 

MuJi) - J2{-uJi)) - {J2M - J2{-^2)) - {J2M - M-^s)) 



9A_ 

2uiUJ2UJo^ V F 
i 



{LJ2U3 -q2-q3){S 



1 



(J2(a;3) - ^2(-C^3)) 



y ^ (a;2a;3 - q2 • ^3) (^S^ + ^ 



2uJ\UJ2Uz 

- {J2{UJ2) + J2{-UJ2)) + WUJS) + J2{-i^3)) 
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^1 = 



L>4 = 



9aY ( a2 1 



9A\^fc,2 1 



^ujiUJ20J^ \ F ) V 2 
{J2M - J2{-uJi)) - {J2M - J2{-0J2)) - {J2M - M-^s))] (A.25) 



Diagrams s4+s5: 



Da 



S {UJ2OJ2, - g2 • ga) 
1 



4 \F J 

J2{(JJ2) - J2{-^2) 



0J\ 



1 



+ 



^2 = 1^ f ^) S'^{U2U2, - q2 ■ qs) 



F J 



J2(a;i) - J2{-(^i) ( 1 



, ,2 



1 ^ 


) 
















) 







^ ^2(^0^2) - J2(-C^2) /^J_ ^ ^\ ^ ^2(0^3) - J2{-UJ^) [l_ ^ l_\ 
OjI \I-jJi UJ3J UJ^ V OUi UJ2) 



J (y) '5^(^2C^3 - q2 ■ qs) 

, .-/2(-'2) + J2(--'2) / 1 



J2M + J2{-^ / J_ _ J_\ 

\a;2 UJ3J 



, ,2 
•^1 



1 1 

0J\ OJ2 



+ 



J2(a;i) + J2(-cui) _ ^ 
Vcc;2 ci;3 



J2(cc;i) - J2(-CJl) _^ J_ 

Lo\ \^2 



+ 













0J2J _ 



3^ /^^y ^2 

^ J2(CU2) + J2(-C^2) ^\ ^ J2(C^3) + J2(-C^3) f]_ ^ ]_ 

\ UJi UJ2 



^2 



3i (gA^ ^2 
'%\f) ^ 



^2(0^1) - J2(-CUl) 1 I 1 



^ J2(C^2) - ^2(-a^2) _ M + ^2(t^3) - J2(-a;3) _ ^ 

a;! \ijJ\ oJz) iol uj2 



(A.26) 



Diagrams s6+s7: 



^1 



^ Qa 



-(a;2 + '^3)A7r 



(a;2 - ujs) {uJs{Jo{uj3) - Joi-ujs)) - uj2{JoM - M-^2))) 
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^2 = 



"8a;i 

_j_9A_ 



= D4 = 
Diagram s8+s9: 



{llI2 - UJ3) i^siM'^a) + M-'^a)) + i^2{M^2) + Jo{-'^2)) 

{Ui + Us) {uJ3{JoM - Jo{-UJ3)) + UJi{Jo{uJi) - Jo{-UJi)) + 3A^) 



- J2(-C^l) / 



+ 



352 



2 



2(g^ + l)^ 

U2UJ3 
1 252^ 



1 



tc;2 



UJ1UJ2 



UJ2UJ3 



J2{cUs) - J2{-U3) f + 1 



OJ3 



UJ2UJ3 



+ 



UJILJ3 ^ 
U1U2 , 



9a\ 



A2 = 7 {-ft) ('^2<^3 - q2 ■ gs) 



FJ 

J2M - J2(-a;i) / 



UJ1UJ2 OJ1OJ3 



UJ2003 



J2M - J2{-u^ /30S2 + 1) ^ 3{S^ + 1) ^ 



UJ2 



UJ1UJ2 



U!lU!3^ 



J2M - J2{-i^3) (S^ + l 



UJ3 



UJ1UJ3 



UJ2'^3 



<U2'-^3 

1 ^ 25^ + 3 ^ 



a;ia;2 



J2(c^l) + J2(-C^l) - 2 J2(0) (S^ ^\ 

I a;ia;2 ^\''^3 ) 



U^2) + J2(--'2) - 2,72(0) + 1 g2^^ + r 



C<;2 V ^I'^S ti^2<^3 

J2(a;3) + J2{-^3) - 2^2(0) + 1 -52 + 1 



OJ3 



0J]0J3 



OU2OU3 



UJ1UJ3 ^ 

C01UJ2 . 



- ( — ) {UJ1OJ3 - qi ■ qs) 



J2M-J2{-^l) fS^ 
J2{^2) - J2{-^2) 



1 s^ + i 252^3^ 



UJ1UJ2 



UJ1UJ3 UJ2^3 



0J2 



CJ1CJ2 UJ2^3 



UJ1UJ3 



J2(a;3) - J2(-C^3) (^{S^ + l) , 3(52 + 1) 



UJ3 



UJ1UJ3 



+ 



UJ2OJ3 ^1^2 I 



(A.27) 
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Da = 



Diagram slO: 



J2M + J2{-ooi) - 2J2(0) / 3 



(jJl \UJ1UJ2 OJiOJs 

J2M + J2i-UJ2) - 2J2(0) fj- 1 

UJ2 \LJ1UJ2 UJ2i^3 i^l^^S 



+ 



^ J2{UJ3) + J2{-UJ3) - 2J2(0) I^J- 



1 2 
+ 



UJ3 



U)iU)3 UJ2OJ3 OJ1UJ2 



2 \ 



J2{UJ2) - J2{-<^2) 1 


+ + 

\UJlUJ2 OJ1UJ2, 0J2^3 

^ 1 1 2 \ 


U2 ' 
J2{0Jz) - J2(-<^3) / 


\UJlUJ2 UJ2UJ3 ijJxUJ^) 

' 1 1 2 \ 



003 



Va;ia;3 0:2003 UJ1UJ2J 



(A.28) 



^4 = 
B3 = 

= 

Diagram sll: 



« 9A 



-{UJ2+UJ3)A^ 

i 9a 



8cui 



{002 - OO3) UJ3{Jo{(^3) + Jo{-(^3)) + 002{Jo{(^2) + Jo{-002)) 



i 9A 



8LO2 



{uji + ujs) {-uj3{Jo{uj3) - Jq{-uj3)) +a;i(Jo(a;i) - Jo{-uJi))) 



(A.29) 



p J (^2^3 - 92 ■ 93) (^'^ + 2 



J2(a;i) - ^2(0) 



CUT 



1 1 

UJ2 <^3 



J2(-a;i) 



^2(-cu2) 1^ 2 



J2(CU2) - J2(0) / 2 ^ 1 

a;2 va;ia;3 UJ1UJ2 

J2M - ^2(0) ^ 2 ^ 1 



<^3 



+2J^(0) 



Miu;2 

1 1 

+ 



<^l<^3 
1 



a;2a;3 
1 

a;2C<^3' 



a;ia;2 OU1OU3 002OU3 



UJi \UJ2UJ3 UJ1UJ2 

J2i-C02) - J2{0) n 
UJ2 

^1 _ ^2(^0^3^2(0) M 
UJ3 \UJi 



1 



UJ1UJ3 
1 



(jJ3 
1 

a;2 
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Va;2 



J2M - ^2(0) / 2 



^3 



+2J^(0) 



1 



+ 



1 



1 \ 



1 



+ 



Va;i UJ2, 
J2i-uJ3) - ^2(0) 



UJn 



1 



0^2 / 



UJ1UJ2 
5 



a;ia;3 a;2i^3 



^4 = - J (y) (c^2C^3 - g2 • qs) + ^) 

J2(c^i) + J2{-u;i) - 2J2(0) _ 1_\ 

\UJ2 UJsJ 



J2M - ■/2(0) 



— ^ — + — 

UJ2 UJ1UJ2 

J2M - J2(0) ^ 2 _^ ^ 



1 \ 



1^21^3 

1 



J 



0^3 

5 



a;ia;2 C(;iu;3 u^2'^3 



+ 



J2(-C^2) - J2(0) _ l_\ 

UJ2 \ ^1 ^3 ^ 

J2(-C^3)-^2(0) 

\a;i a;2 



y j (u;2u;3 - g2 • ^3) (s^ + ^) 
J2M - J2(0) ^ M + ■^2(-c^i) - J2(0) / 



a;2a;3 a;ia;2 coicu^ 



^ J2M - J2{-UJ2) (]_^]_\^ J2M - J2{-0Jz) (]_^]_ 

UJ^ \uJi UJ2 



^2 



-2J^(0) 
_1 



111 
— + + — 



ijJ\ijJ2 



4 VFy 



3^ 
2. 



a;2<^3 



J2(c^l) + J2(-C^l) - 2J2(0) /J_ 

Va;2 

J2(C^2) - J2(0) / 2 ^ 

ci;2 \a;i(x;3 a;ia;2 

1 



1 

<^3 



a;2<^3 

^ ■/2(c^3)-./2(0) / 2 ^ 

CJ3 \a;i6i;2 Ci;ia;3 a;2<^3 



1 ^ J2(-UJ2) - J2(0) (1_ _ 



^2 



M _ J2(-C^3) - ■/2(0) 



^2(0^1) - ^2(0) 

J2M - ■/2(0) / 2 



J_ ^ J_^) ^ J2(-a;i) - J2(0) 



OJ2 OJs 



+ 



0J2 



L\ + - J2(o) / 1 



a;2a;3 UJ1UJ2 <^i(^3 

<^3/ 



1 \ ^ J2(-C^3) - J2(0) _ 



^ J2M - J2{0) 

\U1U2 ' UJ1UJ3 



+ 



UJ2UJ3' 



OJi U2) 
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-2J^(0) 



1 1 

1 



LiJ\iiJ2 OJ\iiJ-i UJ2UJ3 



Diagrams sl2+sl3: add up to 0. 
Diagrams sl4+sl5: 



A2 



.9a 02 J2M - M-^i) 
A4 = S3 = Di = D4 = 



Diagrams sl6+sl7: 



A2 
A, 

A, 



= D4 = 



Diagrams sl8+sl9: 



Ao = 



^3 



= 

A, = 

Diagrams s20+s21: 



1 

2 



5^ + 



1\ CUI+UJ3 



4 V~ ■ 2 
= = 



<^2 



(J2(^i) + J2(-u;i)-2J2(0)) 

(J2('^2) - -^2(-'^2)) 



A, = 



3^ 9a c2 ^2 - Cl^3 

4 Fs'' 

3« 9a c2^1 +'^3 



4 F5 
L)^ = L)4 = 



(J2('^i) - -^2(-a;i)) 



{M^2) - ^2(-<^2)) 



Diagrams s22+s23: 



^1 = -^-^ ['^2(-'"o('^2) - Jo(-'^2)) + u;3{Jo{u;3) - Jo(-a;3)) + 2A, 
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A2 
^4 



^9a_ 
4 F5 

3i 5'A 



^3 = 

Diagrams s24+s25: 



Diagram s26+s27: 



A,, 



B3 
A, 



■ 9a fc2 , i\ M'^^i) - M-'^i) 



"'f^I*^ +2J 

/^C^2 , ■/2(C^2) - J2(-C^2) 

'f^ r + 2 j 

Di = = 



A2 
A4 



4F5l'^ +2 



1 \ a;2 - a;3 



J2(u;i) - J2(-t^i) 



- 2J^(0) 



4F5 +2 



^3 = 

= Di = D^ = 
Diagrams s28+s29: 



^fA + 2J a;f 

^ 9a ( c2 . 1\ ^1+^3 



U}2 



^2(0^1) + ^2(-Wl)-2J2(0) 



J2(a;2) - J2{-UJ2) 



0U2 



-2J^(0) 



[ -^o(^3,C^l) - Kojts, -UJi) ^^^2 _ ^^^2 

Jo(a;i) - Jo(-cc;i) 



{2tl - hM% + 2M^) 



ci;2a;3 - q2 ■ 



-2t3 + 2u;f - M') + /o(t3) 



-2t3 + M2 , A, 



+ 



^2 ~ 1^2 



ii^0(t2, ^2) - i^0(t2^(^2 _ 4^2^^ ^ 4^2^2) 



8a;o 



Jo(a;2) - Jo(-'^2) 



8a;2 
^2^?, - 92 • g2 



(i2-2a;2^) + 7o(t2)if 



to ^2 ~ 1^9 



167r2 



8a;3 



^0(0^3) - ^o(-a^3) . o, .2^ I , ^1 " ^1 
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20 . 



-(Jo(cJi) - M-uJi)){h + 2M^ - 2uf) 



-/o(t3)7^(5ti - 8M\ + ujh - 4MV) 



(QM'' - ts){t3 - cuf) 



(^^2(^3 - q2- qs 



ti - 



8a;3 

(.,-2.I) + /o(.0| + ^ 



■ 9a i (^2 - ^s) 
'F^\8{ts-u;!) 

-{JoM + Jo{-uJi)){ts + 2M^ - 2uf) - Jo(0)2(t3 - 2M^) 

^ {ui2U}s -q2- qz) 

8uJ2{t2 - ^l) 
+ (Jo(6^2) + Jo{-^2)){t2 - 2ujI) + Jo(0)2t2 



-(i^o(^2, ^2) + i^o(t2, -^2)m - 4M% + 4MX) 



{UJ2UJ3 - ^2 ■ 



8cu3{ti - cui) 
+ iMou3) + M-uJ3)){ti - 2cul) + Jo(0)2ti 



^3 = # 



<^i<^3 - gi • g3 

^3 -u;? 



+ ^°^"^^:^°^""^^ t3-2.D+/o(t3) 



8a;i 



8uji 

2 167r2 



{Ko{t2,UJ2) - Ko{t2, -UJ2))t2{t2 - 2M') 

20 . 



8(^2 - ojD 

{JoM - M-0J2)){t2 + - 2ujI) 
2 



3a;2 



--^0(^2)^(5^2 - 8M^i2 + ^^2^2 - 4:M^U2) 
0LU2 



1 



367r^a;2 



(6M^-t2)(t2-^) 



1 



Vsj 16a;2(t2-cj2) 
-(Jo(u;2) + Jo(-^^2))(^2 - 2u;l) - 2Jo{0)t2 



+ 



1 
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-{JoM + M-uJ3)){ti - 2u;l) - 2Jo(0)ii 



Ko{t3,UJi) - Ko{t2, -UI2) 



OU2 



Jq{uJ2) - Jq{-UJ2) 



<JJ2 



{t2 - 2ul) + /o(t2)4t2 



t2 - UJ2 

27r2 



tl — 



<^3 



COS 



Diagrams s30+s31: 



A 



. 9 A / 

1 = «tt;7<^ 



J- / LJ2{Ko{t2,UJ2) - Ko{t2, -UJ2)) 

AF^\(t2-u;'2)'[ 2 
{qi ■ (gi - gs) - 0^10^2) (-4t2M2 + 4M2a;^ + 3i^) 

-(^2 - cJ2)(3t2 - 6^2 + 4:M^uji) 



-{t2-ujl)it2-2M^ + 2uj) 



+/o(^2) 



2a;ia;2i2(2t2 + (^2) + 



3^2 



+A^(t2-a;2)' 



t2 — 

167r2 



912 



'^^^'~''\t2-QM' + l^ul) 



9 



/ a;3(-f^o(^i,t^3) - -f^o(^i, -^^3)) 



(gi • {qi - q2) - uJiUJ3){-AhM^ + AMW3 + 3tl) 

, '^3(^o('^3) - ^o(-'^3)) 



{qi • {qi - ^2) - u;iUJ3){5ti - 2u;l) 
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(A.38) 



-(ti - ujl){ti - 2M^ + 2ujI) 



+Io{ti) 



2uJiUJ3ti{2ti + ujD + 



ti - Lol 



+ gi • (?i - g2)(ll^ic^l - - 2MH\ + ^MHxuI - 2M'^ul + 2t\) 

3^1 



ti - a;| r 2 



167r2 
9 



(ii - 6M^ + 18a;^ 



-(Xo(^2,t^2) - -f^'o(^2, -a;2))t^2^ 



+Io{t2)t2{t2 + Sul) - A^2(t2 - 

1 / ^3(^0(^1,^3) - ^0(^1, -^3)) 

(ti-a;|)n 2 
;-g2 • (gi - 92) + cu2CJ3)(-4tiM2 + 4M^ujI + 3t?) 
-2(ti -cu2)ti(ti - Af2)" 

(-92 • {qi - 92) + uj2UJ3){5ti - 2ujI) 



-2{ti-ujl){h-M^) 



+Io{ti) 



^92 • (qi - 92)(lli?c^3 - 4.tiu;^ - 2MH\ + AMHiuI - 2M^u;l + 2tl) 

Oil 



-2cj2^3^i(2ti + ujI) + 2{ti - ujl){ti - M^){ti + 3ljI) 
4 

-A^{ti - u;3f—q2 ■ {qi - ^2) 



^1 -^3 
167r2 

9a f 



9ti 



3ti 



q2 ■ {qi - 92) (q - 6tiM^ + ITtiWg + GM^uj^) - Auj2UJ^ 



{Ko{t2,UJ2)+ Ko{t2,-UJ2)) 



OJ2 

AF^\2{t2-ujlf 

{ujx + uj3)uj2{-4.t2M'^ + 4.M'^ujI + 3^2) - (^2 - i^l){tl - ^hM'^ + AM'^lo^) 



+ (Jo(a;2) + Jo(-'^2)) -(1^1 + 1^3)^2(5^2 - 2ul) + (^2 - '^2){t2 - 2^^ 



+ Jo{0)t2 G{UJ1 + UJ3)U2 + 2{t2 - U)l) 
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+ {Jo{^3) + M-^s)) (<^i +<^2)<^3(5ti - 2ujI) - {ti - ujl){ti - 2ujI) 
-Jo(0)ti 6{uJi + uj2)ujs + 2{ti - ul 



■ 9 a / 
= I — ^< - 



'a;i(i^o(^3,^^i) - -f^o(^3, -^i)) 



-q3 ■ {q2 + gs) + iOiL03){-4t3M' + AM'ui + 3t|) + 2(^3 - uji)h{h - M') 



Ui{Jo{ui) - Jq{-Ui)) 



+2{h-ujl){h-M^) 



(-?3 • (^2 + ?3) + i^li^z){^h - 2ul) 



i^qs ■ (?2 + q3){ntiu;i - Atscot - 2M% + 4.M%uji - 2M'ujI + 2tl) 

613 



-2uiUJ3h{2h + c^i) - 2(t3 - ooDih - M^){t3 + ?>ujI) 
4 

-A^(i3 - ujif—q2 ■ {q2 + qs) 



167r2 
1 



9t: 



3t 



-qs ■ {q2 + q3){ti - QhM' + llhuji + QM'uji) - AujIuj^ 



{Ko{lJ3) - Ko{-LU3))LLl3tl + {Jq{llI3) - Jo{-UJ3))c03{^tl - 2ujj) 



Di = = 
Diagram s32: 



(A.39) 



Diagram s33: 



B3 

A, 



:9l fg2^^ 



'2F5 



2 

^4 = 



{uji{t3 - cuf) 



{Ko{t3,UJi) - Ko{t3, -Ui)) 



M%{t3 - 2u;j) 



16 



+ (Jo(o;i) - Jo(-a;i))^^ -r^ + h{h)—r^ 



48 



6 

q2 ■ {qi - qs) 



167r2 



18 



(A.40) 
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{Ko{t2,u;2) - Ko{t2, -UJ2)){4 + 2M% - Atiui - S^sM^ + AM'uit2 + SM^) 
+Io{t2Hu;2{tl - 2MH2 + 2uJii2 - AM'ul) + '^^'^^^^ ' ^^2) 



27r2 



qs • {qi - 52) 



{Koiti,LU3) - Ko{ti, -UJ3)){4 + 2MHI - Atjul - 8tiM^ + AM^ulh + SM^cul) 



-{JoM - Jo{-u}3)){ti - 2M\ + 2uiti - AM^ul) 
+Ioih)Au;sitl - 2M\ + 2cvlh - ^MW,) + " ^^'^^^^ " ^^'^ 



(Xo(^3,^^l) - -f^o(^3, -^1)) 



27r2 



U{tl - 2ujft3 - 2MH3 + AW^ujI) 



16 



-{Jo{ui) - Jo(-^i)) 



2>tl - IQM'% + Au'lt^ - Acuf + AM^cuf 
48 



+-io[t3) ■^TT r 



M^{u;i-q2-(qi-q3)) 



6 167r2 18 

- - M-(^2)){t2 + 2ljI) + Io{t2)AuJ2{t2 + 2lJ^) 

UJ2{t2 — UJ2) 



27r2 



16uJ2{t2 - t^l)^ 

+ {Ko{t2,UJ2) - Ko{t2, -UJ2)){-^Ujlt2 + 4- AM^Lul + AM%) + 
- 93 • (gl - g2) 

16a;3(ii-^|)2 

(Xo(^i, ^^3) - Koih, -u;3))(tl + 2MHI - At\ujl - ShM^ + AM^ulh + SM^ul) 

-(Jo(cu3) - M-^3))itl - 2M% + 2ojlti - AM^ul) 

U3{ti - 2M^){ti-ujij 



+Io(ti)Au;3{ti - 2M% + 2ujiti - AM^u^) + 



0J2 — 0J3 



2\2 



{Ko{t3,u;i) + Ko{t3, -u;i))tl{t3 - 2ul) 



-(Jo(cui) + Jo(-cui))-(-4cUi% + ^tl - AM^Lul + AMH3 + Alu^) 
-•^o(0)^(-6a;^t3 + 3tl + AM^ujf - AM%) 

{^l - q2 ■ jqi - g3))(c^i + 
I6(t2-^i)' 

- (^0(^2, c^2) + Ko{t2, -uj2)){AM'^ujlt2 + AMHI - SM^ut - Qtlul + ^2) 
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+ {JoM + Jo{-uj2))-{lQM^u}i - I6MH2 - 4a;^ + 166^2^2 + 3^^) 

+Jo(0)-(-16M2a;^ + I6MH2 - IScu^^s + 3t^) 
3 



UJ1UJ2 - qi- q2 



{Ko{t2, u;2) + Ko(t2, -a;2))3a;|(4M^a;2 ^ ^2 _ 4^2^^^ 



24o;2(t2-o;i)2 
+ (^0(^2) + Jo(-W2))(8M2o;2 - SM^ta + 5a;2% - 2a;^) 
- Jo(0)2(-8M2t2 + SM^uj^ + 3a;^t2) 

(^i-g3- (gi - g2))(c^i +^2) 
16(ti-w|)' 

- (Xo(^i, ^3) + Ko{ti, -oo3))(4^M^oulti + AMHI - m^uj^ - ^IujI + tl) 

o 

2 

+<^o(0)-(-16M=^a;| + IQMHi - 18a;|ti + 3t?) 
o 



+ (Jo(a;3) + Jq{-u}3)){8M^uI - 8MHi + buolti - 2a;^) 
-Jo{0)2{-8M\ + 8M^u}l + 3u}lti)] \ 



—I 



■ aW^i-qi- (92 + gs) 



{Ko{t3,cJi) - Ko{t3, -uJi)){tl + 2MHI - 4:tl(jjl - 8hM^ + AM'^ulh + 8MW^) 
-{JoM - Jo{-uJi)){tl - 2M% + 2ujlh - AM'^ujI) 



+Io{h)Auji{ti - 2M% + 2ujih - AM'uf) + 



uiih - 2M^){t3 - col) 



UJ2{t2 -UjI) 



{Ko{t2,UJ2) - Ko{t2, -UJ2)) 



t2(tl - 2ujlt2 ~ 2MH2 + AM^ul) 



16 



48 



+/o(^2) 



ul{t2 - 2M2) _ ^ a;2(t2 - ^l) , ^2 9t2 - 30M^ + 8 



+ 



M\ul - qs ■ (gi - g2)) 
16a;3(ti-a;i)' 



6 ■ 167r2 18 
(Jo(u;3) - M-uJsWi + 2u;l) + Io{ti)Au;3{ti + 2u;l) 



38 



^4 = -i-^J—-^-^\-{Ko{ts,u,)-Ko{h,-u,)){4MW,+tl-4MH^ 

1.01' [LVi[t3 ~ UJi) 



+ (Jo(a;i) - Jo{-uJi)){t3 - 2ujI) + Io{h)4t3UJi + 



1 



{Ko{t2,uJ2) - Ko{t2, -LU2)){4M^cui + tl- AMH2) 

UJ2{t2 — OjI) 



UJ2{t2 -UjI) 

+ {Jo{uJ2) - M-^2)){t2 - 2ujI) + Io{t2)4:t2UJ2 + 
1 



27r2 



Diagram s34: 



27r2 



(A.41) 



A, 

A2 

B3 
A, 



.9\ 



1 



F'' {q2 + q3- qi? - 
1 



—I 



.9\ 



F^{q2 + q3-qiY-M^ V 2/ ^ ^ 



:9 a 



1 



F'' {q2 + gs - - 



(A.42) 



Diagram s35: Combine with diagrams Rid and 3r. 

4. Rcnormalization diagrams: We now give those contributions from the first order diagrams, 
which come in via renormahzation and thus cannot be obtained from the 1/m expansion 
of the relativistic amphtude. To third order, the renormahzed pion mass, pion Z-factor, 
decay constant and the axial-vector coupling read: 



^ p2 ^4 p2 ' 



Ml 



2M2 A^- 

p2 ^3 + 



F. = 

F^ 



gA\ i 4M2 , gl 



FJo 



F2 



9a 



4F2 



47]-2 



(A.43) 
(A.44) 

(A.45) 

(A.46) 



Here, M, F and gA,o denote the corresponding leading order values in the chiral expansion. 
Moreover, one has to take care of the renormahzation of the nucleon mass when making 
the heavy baryon transformation: 



+ 



3glM' 



V • p V ■ p {v ■ pY 327rF2 
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(A.47) 



The nucleoli Z-factor takes the simple form: 



1 - 



9l 3M^ 
F2 327r2 ■ 



(A.48) 



Note that the familiar In-terms have been treated as described in [^. With that, we 
find for the pertinent diagrams: 
Diagram Rla: 



A, 



Da 



q2 ■ qz] 



3M2 



4 



UJlU}^ ijJ\l02 1-1^2^3 



9\ 



9A 



—I 



9A 

F 

9\ 



{uj2^z - q2- qa) 
8di6(A) A 



2F2 
3M2 



(gl + 2) 



+ 



+ 



UJ1UJ3 UJ1UJ2 UJ2^3 



167r2F2 



9a 



—I 



9a 
F 
1 



1 1 
— + — 

9Ay 

f) 

9\ 



q2 ■ gs) 

+ 



2F2 
3M3 



64:71 

1 / 1 



[UJ1UJ3 - gi ■ gs 
8rfi6(A) 



r2 . 1 





H- ^ ( 


1 








^012 





3M2 / 



+ 



+ 



4 Vco'iCo'3 ci;iCi;2 '-i;2Ci;3 



— z 



167r2F2 c/^ 

F/ 1287r 

1/1 1 
- + — 



rfl6(A) 



UJ1UJ3 \UJi UJ^, 

3 3M2 / 1 



9a 
F 

9l 



UJ1UJ2 



+ 



A^ 
2F2 

1 



{9I + 2) 



1 



V167r2F2 
Diagrams Rlb+Rlc: 



^die{X) A^ 2 
2F^^^^ 



C^2 
1 

0^2 ^^3 

2)1 





1 




002003 V 


C^2 


UJ-iJi 



(A.49) 



^2 
53 

Ai 



. gA 002-003 f 4M^ ^ 

'2F^^;^ 1^ ''^^^ + 32^ + F^ 
.3(7iM3 CU2-CU3 



647rF5 



cut 



gA UJ1+UJ3 /4M2 



2F3 CJ2 

Di = Di = 



+ 



A. 



327r2F2 ' F2 



7 



9l' 



(A.50) 
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Diagram Rid: 



= A^^Di^D^^O (A.51) 
Diagrams Rle+3r+s35: 

■9a 1 f ,,2f'^M\ A \ 



{qi - 92 - qsf - [ \ ^ 2Fy 
^3 = L__,|_M^&+A.^ 



A4 = = = (A.52) 

5. Counterterm amplitude at threshold: 

Finally, we spell out the contribution from the various counterterms to the threshold am- 
plitudes. At threshold, ^/s — m + 2M^, q2 — — (M^, 0) and p2 — (0, 0). The energy 
of the incoming pion and the incoming nucleon can be easily evaluated together with the 
spinor normalization factor A/i. At threshold, the amplitude is fully given by Ai and A2. 
The contribution from the different diagrams reads: 
Diagrams 2a-|-2b: 



. gA M^v ■ pi 
'F^ ul 

^2 = (A.53) 



Ai = -2^^ — — (2C1 + C2 + C3) 



Diagrams Sa+Sb: 



Ai = 

A2 = i^^[2ciM-a;i(c2 + C3)] (A.54) 
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Diagrams 3c+3d: 



^1 
A2 



'-{2Ci + C2 + C3 





Diagrams 3n + So: 



Diagram Bp: 



A2 



Ai = 



Ao 



SM'gA 



_C2_ 

2m 



+ (i26(A) +(i27(A) +(i28(A) 



4M2 
2M2 



9^ - ciio(A) - cii2(A) - die (A) + ^ciis 
Im 1 

_gACA _ ^^^(-^-j _ d^^i^^X) + (ii6(A) - ]-dis 
Im 1 



Diagram 3q: 



Ai = i 



2F^{ui - M) 
.M2(2cui -M) 



—I 



2F3(a;i - M) 



[2di6(A) - dis] 
[2di6(A) - dis] 



Diagram 3r: Combine with Rle. 
Diagram 3s: 



(A.55) 



(A.56) 



(A.57) 



(A.58) 



Ai = 2|^^^[2£iM(2a;i-M) + 24a;i(2M-u;i)+£4M(2M + a;i)] 



.9A M 

^F''LUi - M 



2liUi{2M - ui) + 4(4Mu;i -ul- M^) - ^Mui (A.59) 



Counterterm contribution at threshold from renormahzation diagrams: 
Diagram Rid: 



^1 = 
Ao = 



"^-pF^i6(A) 



(A.60) 



Diagrams Rle+3r: 



Ar 
A2 



.9a 



^9A 



2(a;i - M) 



F^2{u;i-M) 



-4 + 34 , 6 ■ 
9a 

Mes + {2ui-M)h , 2(2u;i-M) , 

H "16(,Aj 



F2 



9a 



{AM) 



42 



Altogether: 



^1 



A,, 



J_i9A 



(2ci + C2 + C3) {-2vpi + '^ 



i?3 



M2 



+ 



m 



m 



+8MV(^^~26(A) + ^27 (A) + ^28 (A)) 

-4M2(dio(A) + rfi2(A)+rfi6(A) 



1 , ^ ^ 



M{2uji - M) + 2£2UJi{2M - ui) 



+4M ( — + cui j + —£3 



4|^(4ciM - 2c.i(c2 + C3) - C4M) 
-2M2(dn(A) + rfi3(A) + -^^18) + 2(Li-M) 



+ ^(4-^4) 



2£iCJi(2M - wi) + £2(4McJi - 00^ - 



(2C2 + C4) 



(A.62) 



(A.63) 
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Figure 1: Tree graphs with insertions from the dimension one La- 
grangian (la-le) and the dimension two Lagrangian as depicted by the 
dot (2a,b). Sohd and dashed hnes denote nucleons and pions, respec- 
tively. 
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Figure 2: Tree graphs of third order in the chiral expansion. The 
filled box denotes an insertion from the dimension three pion-nucleon 
Lagrangian whereas the circle-cross denotes a fourth order mesonic in- 
sertion. 
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Figure 3: One-loop graphs of the tadpole type. 
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Figure 4: One-loop graphs of the self-energy type. 
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Figure 5: Further one-loop graphs of the self-energy type. 
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Figure 6: Fits to tlic total cross sections (solid lines). The dashed and 
dotted lines refer to the second and first order contributions, respectively. 
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Figure 7: Fits to the differential total cross sections for 7r~p 7r+7r~n 
(solid lines) with respect to the kinetic energy and of solid angle the 
outgoing TT"*". T2 = to'2 — refers to the outgoing positively charged 
pion. For further notations, see fig.^. 
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Figure 8: Predictions for the total cross sections up to T^^ = 400 MeV 
(solid lines). For further notations, see fig.^. 
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Figure 9: Predictions for the angular correlation functions at fixed 62- 
For further notations, see fig.|^. 
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Figure 10: Predictions for the angular correlation functions at fixed ip2- 
For further notations, see fig.|^. 
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Figure 11: Predictions for the differential cross sections dajdM'^^ and 
da /dt. For further notations, see fig.||. 
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Figure 12: Chiral expansion. In the upper two panels, the sohd, dashed 
and dotted hnes refer to the contribution up to third, second and first 

order respectively. In the lower panels, the dashed lines show the third or- 
der contribution minus the ones from the dimension three contact terms. 
Similarly, for the dotted lines the loop contributions are in addition sub- 
tracted. 
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Figure 13: Sensitivity to the pion-pion interaction. Considered are 
the two channels Tr+p Tr~^Tr~^n (left) and 7i~p vr^vr^n (right). Upper 
panels: Total cross sections in the threshold region. Middle panels: da/dt 
at = 220 McV. Lower panels: da/dt at = 240 McV. The sohd line 
refers to the standard case, whereas the dashed lines arc obtained by 
setting Is — —70 and keeping all other LECs fixed. 
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